The dust SED in the dwarf galaxy NGC 1569: 
Indications for an altered dust composition? 
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Abstract. We discuss the interpretation of the dust SED from the mid-infrared to the millimeter 
range of NGC 1569. The model developed by Desert et al. f?) including three dust components 
(Polyaromatic Hydrocarbons, Very Small Grains and big grains) can explain the data using a 
realistic interstellar radiation field and adopting an enhanced abundance of VSGs. A simple three- 
temperature model is also able to reproduce the data but requires a very low dust temperature 
which is considered to be unlikely in this low-metallicity starburst galaxy. The high abundance 
of Very Small Grains might be due to large grain destruction in supernova shocks. This possibility 
is supported by ISO data showing that the emission at 14.3 /im, tracing VSGs, is enhanced with 
respect to the emission at 6.7 lim and 850 /xm in regions of high star formation. 



INTRODUCTION 

The dust abundance and properties in galaxies are expected to depend on their metallicity 
and interstellar radiation field (ISRF). In fact, the IRAS colours in dwarf galaxies 
are found to be different from those in spirals |8] implying that the low-metallicity 
environment in these galaxies indeed has an effect on the dust. NGC 1569 is a nearby 
(2.2 Mpc) dwarf galaxy in a post-starburst phase but still exhibiting a high star formation 
rate. It is an excellent candidate to study the dust amount and properties in a low- 
metallicity (12-i-log(0/H)=8.19) galaxy with a high ISRF. 

We observed NGC 1569 at 1200 iim with MAMBO at the IRAM 30m telescope 
and at 450 and 850 jim with SCUBA at the JCMT between 1998 and 2000. A detailed 
description of the observations, data reduction and modelling is given in Lisenfeld et al. 
10] • Here we give a summary of the observations and modelling, compare them to more 
recent results and include a discussion of the spatial variations of the dust SED. 



THE DATA 

We determined total galaxy flux-densities from the maps at 450, 850 and 1200 /Jm 
by integrating them over increasingly larger areas until the cumulative flux-densities 
converged to a final value. However, it was found that low-level emission from NGC 
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FIGURE 1. Left: The interstellar radiation field adopted in the modelling of NGC 1569 (black squares). 
The open triangles are data from Isreal |4], corrected for Galactic foreground extinction which represents 
the bulk of the extinction in NGC 1569. Right: The best fit based on the model of Desert et al. |2]. The 
(blue) crosses are the data used in Lisenfeld et al. ||7|1 and the (red) tiiangles the data of Galliano et al. 



1569 extended over most of the SCUBA field of view of 2 arcmin due to the relatively 
large size of NGC 1569. This low-level emission made a reliable determination of the 
zero-level of the maps impossible. Fortunately, this was not a problem with the larger 
IRAM 1200 jum field. Therefore, we solved the problem by fitting the 850 and 450 
jUm cumulative growthcurve to the scaled growthcurve at 1200 /im, increasing the total 
fluxes at 850 and 450 /im by about 30% in this way. 

The data, together with the IRAS data at 12, 25, 60 and 100 /im taken from the IRAS 
Point Source Catalog, indicate various pecularities in comparison to spiral galaxies. 
The high 25/12 /im ratio indicates that the FAH contribution, dominating the 12 /im 
emission of our Galaxy, must be low with respect to the contribution of Very Small 
Grains dominating at 25 /im. Furthermore, we found a high 25/850 /im ratio and a 
relatively flat spectrum in the (sub)millimeter range. 



MODELLING THE DUST SED 

We used the model of Desert et al. 10] (hereafter DBP90) in order to fit the data. This 
model includes three components: (i) big grains that are in equilibrium with the ISRF, (ii) 
stochastically heated Very Small Grains (VSGs) and (iii) Polyaromatic Hydrocarbons 
(PAHs). Since we only have a limited number of the data points we did not try to exploit 
the full parameter range of this model and only changed the relative contribution of 
the three components. We used the ISRF shown in Fig. 1 (left) derived from ultraviolet 
to infrared observations of NGC 1569 and corrected for Galactic extinction which is 
globally much more important than intrinsic extinction in NGC 1569. 

In Fig. 1 (right) we show the best fit result, together with our (sub)millimeter data 
and the IRAS data points. As stated above, the low 12/25 /im ratio indicates a relatively 



low contribution of PAHs with respect to the VSGs. Within the limited data, we are 
not able to constrain the contribution of PAHs any further and therefore do not include 
them in the fit here. The contribution of VSGs with respect to big grains is enhanced 
by a factor of 5.3 in this fit with respect to the Galactic value. Varying the ISRF in a 
reasonable way, we found that enhancements of the VSGs between factors of 2 and 7 
give good fits to the data. In our fit the VSGs are responsible both for the mid-infrared 
and for the (sub)millimeter emission. Their broad emission spectrum is due to (i) the 
broad temperature range and (ii) the fact that the frequency dependence of their dust 
emissivity is assumed to be (DBP90). 

In this figure we also include data by Galliano et al. [6]. The agreement is generally 
good with the exception of the data point at 60 /im and the SCUBA points which are 
lower than ours. The difference between the two sets of SCUBA data can be explained 
by the aperture correction that we applied. 

As an alternative we tested a multi-temperature model. Because of the limited amount 
of data points we were only able to constrain three temperature components. In spite of 
the simplicity of the model, the results are expected to give a reasonable upper limit to 
the cold dust temperature. We needed very cold dust temperatures in order to explain the 
relatively high flux and the flatness of the spectrum in the (sub)mm range. The best fit 
to the data was achieved with a cold dust temperature of 7 K. A dust temperature of 1 1 
K still yielded acceptable agreement with the data within the error limits. Such low dust 
temperatures can only be achieved in clouds with a very high opacity and a low filling 
factor [6]. However it needs to be tested whether the overall emission of such clouds in 
the high ISRF of NGC 1569 is in agreement with the observed dust SED of NGC 1569. 



SPATIAL VARIATIONS OF THE DUST SED 

Haas et al. [3] compared ISO images at 6.7/im and 14.3/im to SCUBA images at 850 
jum. They found that the emission at 14.3;Um in NGC 1569 is enhanced with respect to 
6.7/im and 850/im in regions of high dust emission, corresponding to regions of high 
star formation. A similar enhancement was found in the nuclear region of NGC 6946 
and in the starbursting overlap region in Arp 244. In spiral galaxies like NGC 891 or 
NGC 7331, on the other hand, there is a very close correlation between the emissions at 
6.7/im, 14.3/im and 850/im. 

This indicates that there is a relation between the star formation intensity and the 
abundance of VSGs which are dominating the emission at 14.3/im. A possible interpre- 
tation is that the supernova shocks produced by the star formation destroy the large dust 
grains by grain-grain collisions |1, 5]. 

The comparison of the maps at different wavelengths in f3\ show that for NGC 1569 
the ratio between 6.7/im and 850/im is constant with a value of 1 .08 ±0.30 (and the ratio 
between 7.7/im, coinciding with a peak in the PAH emission, and 850/im is similarly 
constant with a ratio of 1.80 ±0.60) whereas the ratio between 14.3/im and 6.7/im is 1 
in low emission regions and up to 10 in high-emission regions. In Fig. 2 we fit the two 
sets of data points with the model of DBP90. The emission in low-emission region is 
very similar to the emission of our Galaxy and can be fitted by a low ISRF and solar 
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FIGURE 2. Left: Fit of the model of DBP90 for an ISRF 20 times lower than the field in NGC 1569 
and for relative abundances of PAHs, VSGs and big grains as in the solar neighborhood. The crosses refer 
to the data for low emission regions in NGC 1569 given in Haas et al. |3]. Right: Fit of the model of 
DBP90 for the ISRF as in NGC 1569 (Fig. 1, left). The relative abundances in comparison to the solar 
neighborhood are 3 times lower for PAHs and 6.5 times higher for VSGs. The crosses refer to data for 
high emission regions in NGC 1569 given in fs"] , where we have adopted a relative enhancement of the 
14.3/im emission of a factor of 5. The absolute flux scales are arbitrary. 



neighborhood abundances of VSGs and PAHs. In regions where the emission at 14.3/im 
is increased with respect to 6.7jUm, V.VjUm and 850 /im, an ISRF as in NGC 1569, an 
enhanced abundance of VSGs and a decreased abundance of PAHs are able to explain 
the observations. 
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